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Abstract: 

Vector-boson pair production ranks among the most important Standard-Model benchmark 
processes at the LHC, not only in view of on-going Higgs analyses. These processes may also help 
to gain a deeper understanding of the electroweak interaction in general, and to test the validity 
of the Standard Model at highest energies. In this work, the first calculation of the full one-loop 
electroweak corrections to on-shell W-boson pair production at hadron colliders is presented. We 
discuss the impact of the corrections on the total cross section as well as on relevant differential 
distributions. We observe that corrections due to photon-induced channels can be amazingly 
large at energies accessible at the LHC, while radiation of additional massive vector bosons does 
not influence the results significantly. 
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1 Introduction 



With the recent start of LHC operation hard scattering processes became accessible with par- 
ton energies up to several TeV and, indeed, jet-jet invariant masses up to 5 TeV have been 
observed [TJ[2] in the first round of data taking. With the significant increase of the integrated 
luminosity and the doubling of the beam energy anticipated in the next few years, electroweak 
processes such as the production of lepton or gauge-boson pairs with invariant masses of several 
TeV will become accessible in the near future. These reactions may on the one hand allow for 
precise tests of the Standard Model (SM) and on the other hand potential deviations may point 
to "physics beyond the SM". Indeed, the observation of anomalous couplings of quarks, leptons 
or gauge bosons might well be the first signal of "New Physics". As an alternative one may 
look for peaks or distortions in the W, Z or fermion spectra resulting from the decays of new 
massive particles. Clearly, any well-founded claim of "physics beyond the SM" must be based 
on precise measurements combined with similarly precise calculations, possibly at the level of 
several percent. 

Obviously this requires theory predictions which must include next-to-leading order (NLO), 
possibly even next-to-next-to-leading order (NNLO) perturbative QCD corrections. Considering 
the smallness of the weak coupling, a w , the need for electroweak (EW) corrections is less obvious. 
In a first step one might try to absorb the dominant corrections arising from the running of the 
fine-structure constant from low scales to Mz and from the p-parameter [3] in properly chosen 
effective couplings. However, this approach is only justified for energies of order J% or Mz- 
As is well known from earlier investigations [4l417j. EW corrections increase with the squared 
logarithm of the energy, and may reach several tens of percent for energies accessible at the 
LHC. In view of their strong dependence both on energy and scattering angle, they will induce 
significant distortions in transverse-momentum and rapidity distributions and consequently may 
well mimic "New Physics" and hence must be carefully taken into account. 

Most of the investigations along these lines have concentrated on the issue of Sudakov 
logarithms, i.e. terms that are at n loops enhanced proportional to a™ log 2n ~ m (s/M£t) with 
m = 0, 1, . . . , 2n, corresponding to leading, next-to-leading etc. logarithmic (N m LL) enhance- 
ment. This line of research was motivated by the necessity of including at least the dominant 
two-loop terms, once one-loop corrections exceed the 20 — 30% level. Using evolution equa- 
tions, originally derived in the context of QED [18H21] and QCD, four-fermion processes have 
been studied up to N LL [T lfTIHli] . W-pair production in electron-positron and quark-antiquark 
annihilation up to N 3 LL [22U23] . Employing diagrammatic methods or the framework of soft- 
collinear effective field theory most of these results were confirmed in the NLL and NNLL ap- 
proximation proMS]. 

As stated above, the intermediate energy region, up to approximately 1 TeV, will be explored 
with high statistics, and the complete one-loop corrections may become relevant. This includes 
real and virtual photon radiation and terms suppressed by M^/s, where s denotes the partonic 
center-of-mass (CM) energy squared. Up to now these complete one-loop corrections have been 
evaluated for tt [25H29] (for earlier studies see [HQlEI]), bb [32], as well as for gauge-boson plus 
jet production, the latter in the on-shell approximation [33 -38J and including W- or Z-decays to 
leptons with all off-shell effects consistently taken into account [39|H0]. 

In contrast, the full one-loop EW corrections to vector-boson pair production at the LHC are 
not yet available, despite their great phenomenological interest. These processes contribute the 
most important irreducible background to the production of a SM Higgs boson in the intermediate 
mass range. A profound theoretical understanding of the underlying physics will also allow 
a precise analysis of the non-abelian structure of the EW sector, in particular of the vector 
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boson self interactions. At the LHC these will be explored with high precision, exploiting the 
high luminosity in combination with the highest possible CM energies. Moreover, vector-boson 
pair production processes are well suited for putting experimental constraints on the existence 
of anomalous trilinear and quartic gauge couplings, since their phenomenological effects are 
expected to be sizable at large invariant masses of the vector-boson pairs accessible at the LHC. 
Considering searches for supersymmetry (SUSY) at the LHC, all vector-boson pair production 
channels constitute important backgrounds to signals with leptons and missing transverse energy, 
as e.g. predicted for chargino-neutralino pair production. 

During recent years, a great effort has been made to push the accuracy of the theoretical 
predictions to a new level. The NLO QCD corrections, including the leptonic decays of the 
vector bosons, have been studied by several authors and are implemented in Monte Carlo pro- 
grams |41H49] . The results have been matched with parton showers and combined with soft 
gluon resummations [50] to improve the predictions for vector bosons produced at small trans- 
verse momenta. The loop-induced channels gg —> V1V2, formally a second order effect in QCD, 
nevertheless play an important role at the LHC due to the enhancement of the gluon luminosity 
and have also been studied extensively |51H54j . Specifically, concerning background estimates to 
Higgs production at the LHC, the contribution of the gluon- induced channel to W-boson pair 
production amounts to 30% after experimental cuts [54]. 

Electroweak corrections for gauge-boson pair production at hadron colliders, however, have 
only been evaluated in the high energy limit |55H57j . This is in contrast to the reaction 
e + e~ — > W~W + , where the complete one-loop calculation has long been available [58H5T] . and, 
following the demands of LEP experiments, the full one-loop correction for electron-positron 
annihilation into four fermions (including resonant and non-resonant amplitudes) has been cal- 
culated [52HE5]. 

To arrive at predictions which satisfy the needs of the next round LHC experiments, we em- 
bark on the full one-loop corrections for W-pair production in proton-proton collisions, where the 
gauge bosons are treated as stable particles. For completeness we also recalculate the well-known 
NLO QCD corrections and, furthermore, discuss two competing processes: W-pair production 
through 77 collisions, and through the quark-loop-induced gluon-fusion process. Employing the 
MRST2004QED PDF set [66] for the 77 luminosity, we observe a surprisingly large contribution 
at large invariant W-pair mass and a pronounced peaking at small scattering angles. In fact, 
for small angles and high energies this parametrically suppressed reaction is comparable to the 
qq induced reaction. Even for large angles (i.e. large px) it still amounts to order of 5% and 
is therefore comparable with the genuine EW corrections. Gluon fusion, in contrast, exhibits a 
fairly smooth angular distribution and is typically of order 5 to 10 percent in the whole kinematic 
region of interest. Obviously all these effects must be taken into consideration. 

The layout of this paper is as follows: In Section 2 we start with a qualitative discussion of 
W-pair production at leading order (LO). We introduce the kinematic variables and recall the 
dominant features of W production: the dependence on the W transverse momentum, W-pair 
invariant mass and rapidity. These results will illustrate the large kinematic region accessible 
at the LHC and the marked differences between the qq-, gg- and 77- induced channels, as far 
as transverse-momentum and angular distributions are concerned. Section 3 will be devoted 
to a detailed discussion of one-loop EW corrections to the reaction qq — > W~W + . We spec- 
ify the renormalization scheme, our input parameters, the treatment of QED corrections with 
the associated soft and collinear singularities, their absorption in the PDFs, and the numerical 
evaluation of hard photon radiation. As stated above the rate of W-pair production from 77 
collisions is surprisingly large and the resulting rates and distributions are studied in Section [H 
The cross sections for qq — > W~W + including EW and QCD corrections are convoluted with 
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PDFs, and predictions for various distributions are shown and compared with those from 77 and 
gluon- fusion processes. For future tests and cross checks tables are presented which precisely 
specify the results at various kinematic points. In Section [5j we discuss the phenomenological 
impact of additional massive vector bosons in the final state. Section [6] contains our conclusions. 

2 W-pair production: the leading order 

At lowest order, C(a 2 ), W-boson pair production at the LHC is dominated by quark-antiquark 
annihilation, 

m 3 -> w-w+ , (2.1) 

where the qk denote one of the light quark flavours, qu = u, d, s, c and b. The corresponding 
diagrams are shown in Fig. Q] for the uu channel. Since there is a finite probability of finding a 
photon in a proton we in addition take into account the photon-induced tree-level contributions 

77 -»■ W~W+ (2.2) 

(see Fig. [2]). In the numerical analysis the MRST2004QED PDF set |66j is used for the photon 
density. Although suppressed by two photon-PDF factors, this channel may give rise to poten- 
tially large contributions at high invariant masses. For completeness, we also re-calculate the 
effects due to the loop-induced channel 

gg -> W~W+ , (2.3) 

which adds a relative correction of about 10% to the LO qq cross section [52J. 

Due to the large CM energies accessible at the Tevatron and the LHC, we can safely neglect 
the lepton and light quark masses (i.e. all but the top-quark mass) unless they are needed to 
regularize infrared (IR) singularities in the real and virtual contributions, as will be discussed 
later. However, the bb channel receives contributions from top-quark exchange in the t channel 
that for consistency have to be taken into account in the Born approximation and in the NLO 
calculation. In this work we assume a block-diagonal CKM matrix with T4b = 1- Therefore, the 
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cross sections (|2.ip are proportional to Yl"k=i ^ifc^fcj = ^ or = 1>^> anc ^ ^ * s sufficient to 
only consider flavour-diagonal processes. Consequently, at the parton level the complete one-loop 
calculation has to be carried out independently for the following three different channels, 

uu^W~W+, (2.4a) 
dd^W~W + , (2.4b) 
bb^W~W + . (2.4c) 

Let us, in a first step, evaluate the LO prediction for qq-, 77- and gg-induced processes. The 
cross sections for W-pairs as a function of pS 1 , corresponding to a cut on the transverse momenta 
of the W bosons, and as a function of M^y, corresponding to a cut on the mass of the W-pair, 
are shown in Figs. [3] and HI respectively. The input values and setup we use for our computation 
will be specified in Section 14.11 together with the rapidity and pp cuts employed throughout. 
Given an integrated luminosity of 20 fb _1 at 8 TeV and 200 fb _1 at 14 TeV, transverse momenta 
(invariant masses) up to 0.35 TeV (1 TeV) can be explored at 8 TeV and up to 0.75 TeV (2 TeV) 
at 14 TeV. Here, we assume that events with one W boson decaying into ev e or \iVn and one 
decaying hadronically can be detected with 30% efficiency and require more than 100 detected 
events. 

When comparing the p™* and dependence of qq-, 77- and gg-induced processes, marked 

differences can be observed. For large the rates of 77- and gg-induced reactions are quite 

comparable. However, the cut on reduces the relative effect of the photon-induced channels 
significantly. This is shown on the right hand sides of the respective figures, where the cross 
sections of the gg- and 77-induced reactions relative to the gg-induced process (denoted as 5 gg 
and <5 77 , respectively) are plotted. Let us first concentrate on the LHC at 14 TeV (LHC14). 
The relative rate for W-pairs from gluon fusion amounts to less than 10% for small pt and 
decreases with increasing j?t (Fig. [3]). The parametric suppression of gg — > W~W + by a 2 is 
compensated by the large gluon luminosity for small r = s/s which, however, dies out with 
increasing r. A qualitatively similar behaviour is observed for a cut on Mww (Fig- El)- The 
relative importance of the photon-induced process, in contrast, increases with increasing p™*, 
and even more so with M^y. This behaviour can be traced to the energy dependence of the cross 
section for the partonic subprocess 77 — > W~W + , which approaches a constant value in the high- 
energy limit, <t 77 ~ 8Tra 2 /Myj, with an extremely strong peaking in the forward and backward 
directions. The cut on p^, when compared to the one on Mww thus leads to a significantly 
stronger reduction of the event rate. Nevertheless, a 10% admixture of photon-induced W-pairs 
remains for p^ above 800 GeV. Qualitatively similar statements are applicable to the LHC at 
8 TeV (LHC8). However, the relative contributions of both gg- and 77-induced reactions are 
significantly smaller, and become completely negligible for the Tevatron. In the right three plots 
of Figs. [3] and H] we also anticipate some of the results on one-loop QCD (5qcd) and EW (5ew) 
corrections. For transverse momenta around 800 GeV, accessible at LHC 14, large negative EW 
corrections amount to nearly —50% and are the main subject of this paper. In fact, they are 
quite comparable to the positive QCD corrections, which are displayed in the same figures. Note 
that the latter strongly depend on the cuts on real jet radiation, as will be discussed at the end 
of Section 13.31 

As indicated above, the difference between the p^ ut and dependence of the cross sections 

can be traced to the different angular distributions of the W bosons in the W-pair rest frame, 
which for 77 is strongly peaked in the forward and backward directions. The marked differences in 
the angular distributions are illustrated in Fig. [5l where the distributions in the rapidity difference 
Ai/ww = 2/w~ ~~ Vw+ (which, for fixed M^w, corresponds to the angular distribution in the W- 
pair rest frame) are shown for the mass intervals [2Af\y, 00] and [1000 GeV, 00], respectively. 
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In the first case, the cross section is dominated by central events with W-pairs of low invariant 
mass and relatively isotropic angular distributions. All three components are peaked at small 
Ayww- Considering their ratios, displayed on the r.h.s. one observes the stronger preference of 
the gg process for small Ayww and, conversely, the enhancement of the 77 process for large 
Ayww- Considering events with large invariant mass, M\vw > 1 TeV, the distribution of the 
three components exhibits a drastically different behaviour. The rate of the gluon-fusion process 
is fairly small and its angular distribution does not exhibit any pronounced structure, resulting 
in a flat rapidity distribution. As expected, W-pair production is dominated by qq annihilation 
with its strong peaking at small t, i.e. in the forward and backward directions. This is reflected 
in the peaks of Ayvvw around ±2.5. The contribution from 77 fusion remains sizable, with 
a very pronounced peaking for large rapidity difference, corresponding to a highly anisotropic 
angular distribution in the WW rest frame. The EW corrections exhibit a completely different 
behaviour. For large Ay\yw> corresponding to small t they are negative but less than 10%. In 
the Sudakov regime, however, corresponding to small Ayww and thus large t, the Sudakov effect 
comes into play and negative corrections of more than 30% are observed. In total, the sum of 
these corrections, denoted £<5 in the right displays of Fig. [5j will lead to a dramatic distortion 
of dcr/dAyww with corrections varying between +30% and —30% for large and small |Ayww|> 
respectively. Note that such a behaviour could well be misinterpreted as a signal for anomalous 
couplingso A less pronounced, but qualitatively similar behaviour is observed for LHC8 with 
the mass intervals [2Mw,oo] and [500 GeV,oo] (Fig. [6]). 



3 Radiative corrections 

In addition to the tree- level contributions, a full 0(a 3 ) analysis of the W-boson pair production 
cross section requires the inclusion of the one-loop EW corrections as well as the contributions 
due to radiation of one additional bremsstrahlung photon. Thus, the corresponding total partonic 
cross section at NLO may be written as 

-ij->WW(7) _ -ij-s-WW , 4j-+WW , -ij->WW7 /q -. \ 

a NLO — a LO +°"virt +°"LO ' \ 6 - L ) 

with i, j = u, d, s, c, b, 7 and the corresponding antiparticles, and the different partonic contri- 
butions are given by 



>WW7 



wh I d<I> ( w_w+ )EEE0- M oi 2 + 2Re (- M s-^i)] , (3-2) 

col spin pol 

J dd>(W-W+ 7 ) EEE^I 2 - (3-3) 



V y 

1 1 



col spin pol 



Here, A^o, Mi and A4 7 denote the corresponding Feynman amplitudes of the tree- level, one-loop 
and bremsstrahlung contributions, respectively, j d^ (final state) is the Lorentz invariant phase- 
space measure of the final-state particles and s is the partonic CM energy. The normalization 
factors are given by N q q = 36 and A^ 77 = 4. Defining the momentum fractions x a and X), of the 
initial-state hadrons carried by the incoming partons, the hadronic NLO cross sections at the 



The interplay of logarithmic EW corrections and anomalous trilinear gauge-boson couplings has been studied 
in Ref. [67] for WZ and WW production. 
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Figure 3: Left: Total cross sections as a function of the cut on the W-boson transverse momentum 
at the LHC14 (top), LHC8 (center) and the Tevatron (bottom). The corresponding relative 
rates w.r.t. to the qq channel are shown on the r.h.s., together with the rates of QCD and EW 
corrections. See text for details. 
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Figure 4: Left: Total cross sections as a function of the cut on the WW-invariant mass at the 
LHC14 (top), LHC8 (center) and the Tevatron (bottom). The corresponding relative rates w.r.t. 
to the qq channel are shown on the r.h.s., together with the rates of QCD and EW corrections. 
See text for details. 
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Figure 5: Differential LO cross sections for the W-boson rapidity gap with default cuts (top) 
and with a minimal invariant mass of 1000 GeV (bottom) at the LHC14. On the right-hand- 
side, the corresponding relative rates due to photon- and gluon-induced channels w.r.t. the qq- 
contributions are shown, as well as the EW corrections. See text for details. 



LHC is given by a convolution of the partonic cross section with the parton distribution functions 
(PDFs), 



r NLO WW(7) = / dT f — ^fi/p( X a,^)fj/p(Xb,^)^L(T WM ( T ^ ! J 'i ] ■ (;>> - l ) 

Jrn Jt x b 4 . 



where the hadronic CM energy s is related to s via s = ts, with r = x a Xb. The kinematic 
production threshold of a W-boson pair in the final state is reflected in the choice of the lower 
integration boundary tq = 4M^/s, corresponding to a minimal partonic CM energy of so = T o s > 
and /Up denotes the factorization scale. 

We have performed two completely independent calculations of the EW corrections, and the 
results displayed in Tables [THU have been reproduced by both of them. 
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Figure 6: Differential LO cross sections for the W-boson rapidity gap with default cuts (top) 
and with a minimal invariant mass of 500 GeV (bottom) at the LHC8. On the right-hand- 
side, the corresponding relative rates due to photon- and gluon-induced channels w.r.t. the qq- 
contributions are shown, as well as the EW corrections. See text for details. 



3.1 Virtual corrections and renormalization 

The virtual one-loop corrections receive contributions from self-energies, triangles and box dia- 
grams. In the first approach the diagrams are automatically generated with FeynArts 3 . 5 [68|169j 
and FormCalc 6 . 1 |70p71| is used to calculate and algebraically simplify the corresponding am- 
plitudes. Afterwards, the multiplication with the Born-level amplitude, as well as the summation 
(averaging) over polarisations, spins and colours is performed completely analytically within the 
FormCalc framework. 

In an alternative approach we used the program QGraf [72] to generate Feynman diagrams. 
The implementation of Feynman rules, the reconstruction of the Dirac structure (as well as 
the evaluation of Dirac traces) and the calculation of squared matrix elements is carried out 
analytically using the computer algebra program FORM [73]. To avoid numerical instabilities, 
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potentially small Gram-determinants, which occur in the tensor reduction of particular four- 
point-functions are cancelled at the analytical level using FORM. 

We use two different analytical implementations of the Passarino-Veltman algorithm [71] 
(based on Mathematica and FORM, respectively) to reduce tensor coefficients to scalar inte- 
grals. The reduction was also tested against the numerical approach implemented in the Loop- 
Tools 2.5 |70p75] library, which is used for the evaluation of the scalar one-loop integrals. 

The ultraviolet (UV) divergencies that arise in the computation of the one-loop diagrams are 
treated in dimensional regularization going from 4 to D = 4 — 2e space-time dimensions, where 
the UV divergencies appear as single poles in the small complex parameter e. After adding the 
counterterms in a proper renormalization procedure, the poles vanish, and the limit e — > can 
be taken to obtain physical results. 

Our results are based on the on-shell renormalization scheme defined in the following. 



3.1.1 On-shell scheme 

Our choice for the renormalization prescription is the on-shell (OS) renormalization scheme 
as specified in Ref. [76]. Instead of defining the electromagnetic coupling constant a in the 



Thomson-limit, however, we work in the scheme where a is derived from the Fermi-constant 
via 

^. = ^^(1-4). (3-5) 

In this scheme, the weak corrections to muon decay Ar are included in the charge renormalization 
constant 5Z e by the replacement 

Ar 

5Z e \ — > 5Z e \ = 5Z e \ , n . (3.6) 

Hcr(O) la G M L lo(0) 2 

in the calculation of the counterterm contributions (see, e.g., Ref. |77J). As a consequence, 
the EW corrections are independent of logarithms of the light-quark masses. Moreover, this 
definition effectively resums the contributions associated with the running of a from zero to the 
weak scale and absorbs some leading universal corrections oc G^m^ from the p parameter into 
the LO amplitude. 



3.2 Real corrections 

In a second step, the diagrams contributing to the bremsstrahlung amplitude .M 7 (see Fig. [7]) 
must be considered. In the first approach the amplitudes are generated with FeynArts and ana- 
lytically squared within Mathematica using FeynCalc [78] , We use MadGraph 70] and FormCalc 
for internal checks, however, the computational performance of the FeynArts/FeynCalc-based 
code turns out to be more efficient. Alternatively, the corresponding Feynman-diagrams are 
generated with QGraf and the corresponding amplitudes, as well as the squared matrix ele- 
ments, are evaluated analytically with FORM. In both approaches, the numerical evaluation of the 
bremsstrahlung contributions is carried out in FORTRAN using the VEGAS [80j algorithm. 

In the computation of the real corrections, care has to be taken since the phase-space integral 
over |.M 7 | 2 exhibits IR singularities in phase-space regions where the photon is radiated collinear 
to an initial-state quark or becomes soft, i.e. the energy of the photon goes to zero. We apply 
the well-known technique of the two-cut-off phase-space slicing to analytically carry out the 
phase-space integration over the soft and collinear singularities by using small mass regulators 
m q and A for the light quarks and the photon, respectively, and exploiting universal factorization 



10 



Figure 7: Generic bremsstrahlung diagrams for the process qq — > W W + 7- 



properties for the squared bremsstrahlung amplitudes in the soft and collinear limit (see, e.g., 
Ref. [81]). Accordingly, the IR singularities appear as lnm 9 and In A terms in this particular 
mass regularisation scheme. Note that the hierarchy A <C m q has to be respected carefully in the 
evaluation of the IR-singular terms. To separate the soft and collinear phase-space configurations 
from the hard bremsstrahlung, one imposes a cut on the photon energy, E^ < AE -C vl and 
a cut on the angle between incoming (anti-)quark and photon, 9 qi < A9 <C 1. The residual 
phase-space integration corresponding to hard, non-collinear photon emission, i.e. E~ > AE and 
9 qi > A9, can safely be carried out without IR regulators, allowing for an efficient numerical 
evaluation using Vegas. Adding the soft, collinear and hard contributions, the dependence on 
the slicing parameters AE and A9 cancels out in the computation of IR-safe observablesJl 

According to the Bloch-Nordsiek theorem [82 j , the soft singularities emerging from real radi- 
ation cancel against corresponding contributions from the virtual corrections related to photon 
exchange between on-shell legs in loop diagrams. However, the initial-state collinear singularities 
survive and have to be absorbed in the renormalized PDFs, where we apply the MS factorization 
scheme as described in Ref. [83]. Since we do not include the leptonic decays of the final-state 
W bosons in this work, we do not have to deal with subtleties concerning infrared safety arising 
from the event-selection procedure for photon radiation off charged leptons in the final state. 

2 Technically, the dimensionless quantity S s = 2A_B/\/I is used as slicing parameter instead of AE. 
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3.3 NLO QCD contributions 



In addition to the EW corrections, we have also recalculated the NLO QCD corrections to 
W-boson pair production, using the same setup as detailed above. The corresponding real 
corrections exhibit additional partonic channels with one gluon in the initial state, namely 

qg ^ W~W + q, (3.7a) 
gq W~W + q, (3.7b) 

where the contributions gb ->■ W^t* ->■ W~W + b and gb ->■ W + t* — > W + W~ b formally 
contribute to associated production of a W-boson and a potentially resonant top-quark at leading 
order and thus have to be treated separately. We exclude those channels from our analysis by 
discarding events with a b-quark in the final state, assuming a 100% tagging efficiency. 

In the numerical analysis one finds that the px distribution of the W-boson is plagued by huge 
NLO QCD X-factors of a few hundreds of percent, resulting in potentially large uncertainties in 
the theory predictions. To circumvent this problem, we follow the strategy proposed in Ref. |39j 
and veto events with a hard jet recoiling against a hard W, because these signatures belong 
to W+jet production at LO rather than being a correction to W-pair production. Technically, 
we discard events where the transverse momentum of a visible jet (with px jet > 15 GeV and 
|2/jet | < 2.5) is larger than half of the highest W-boson px- (The important issue of giant QCD 
K-i actors has also been discussed in Ref. [83].) 

4 Numerical results 

In this section we present numerical results for total cross sections and differential distributions for 
the process pp —> W~W + +X at the Fermilab Tevatron for a total CM energy of y/s = 1.96 TeV, 
as well as for the process pp -> W" W+ + X at the CERN LHC with CM energies of y/s = 8 TeV 
(LHC8) and y/s = 14 TeV (LHC14), respectively. We discuss in detail the phenomenological 
implications of the EW corrections. In the following the relative corrections 5 are defined through 
cnlo = (l + <5) xolo- The photon- and gluon-induced channels (|2.2j) and ([2.3)1 . although formally 
contributing at leading order, are also considered as relative corrections <5 77 and <5 gg , respectively. 

4.1 Input parameters and event selection 

We use the following SM input parameters for the numerical analysis, 

G M = 1.16637 x 10~ 5 GeV -2 , 

M w = 80.398 GeV, M z = 91.1876 GeV, 

M H =125 GeV, m t =173.4 GeV, (4.1) 

which are taken from Ref. [85]. In the on-shell scheme applied in our computation, the weak 
mixing angle cos 2 9 W = M^/M^ is a derived quantity. For the computation of the process 
qq -)• W~W + and its EW radiative corrections, we use the MSTW2008LO PDF set [86] in 
the LHAPDF setup [87] . In order to consistently include 0(a) corrections, in particular real 
radiation with the resulting collinear singularities, PDFs should take in principle these QED 
effects into account. Such a PDF analysis has been performed in Ref. [66J, and the 0(a) effects 
are known to be small, as far as their effect on the quark distribution is concerned [88]. In 
addition, the currently available PDFs incorporating 0(a) corrections [66] include QCD effects 
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at NLO, whereas our EW analysis is LO with respect to perturbative QCD only. For these 
reasons, the MSTW2008LO set is used as our default choice for the qq process. 

W-pair production through gluon fusion is considered at leading order, hence the same 
leading-order set [86J is implemented. For W-pair production through the partonic subpro- 
cess 77 W~W+ we use the MRST2004QED set [66]. Note that this set has been also used in 
Ref. [89] to predict the contamination of the Drell-Yan process by the 77 — > contribution. 

For the computation of the QCD corrections, however, the MSTW2008NLO set is adopted, 
corresponding to a value of a s (Mz) = 0.1202. As explained above, the CKM matrix can be 
assumed to be diagonal. The renormalization and factorization scales are always identified, our 
default scale choice being the phase-space dependent average of the W-boson transverse masses 

Mr = = mr = g {\J M w +Pt,w- + \/ M w +#r,w+) ( 42 ) 

for the evaluation of the QCD as well as the EW corrections^ A similar scale choice was taken 
in Ref. [SB] for the computation of the EW corrections to four-lepton production at the LHC. 
In our default setup, we require a minimum transverse momentum and a maximum rapidity for 
the final-state W bosons, 

p T>w ± > 15 GeV, |y w ±| < 2.5, (4.3) 

to exclude events where the bosons are emitted collinear to the initial-state partons. As far as 
QCD corrections are concerned, the jet veto introduced before is only applied to visible jets with 

PTjet > 15 GeV , M<2.5 (4.4) 

to ensure infrared safety. 

4.2 Integrated cross sections 

We start the discussion of our numerical results with the presentation of integrated cross sec- 
tions at the LHC and Tevatron, evaluated with our default setup. Table Q] shows integrated LO 
cross sections, together with the (relative) EW corrections, contributions from 77, gluon fusion 
and QCD corrections. Moreover, relative corrections due to massive-boson radiation (denoted 
as 5wwv) are included, which will be discussed in the following chapter. In addition, the cor- 
responding results for the total cross sections without any event-selection cuts and without a 
dedicated jet veto are shown in Table [2] as a benchmark. While our default cuts reduce the LO 
(/(^-induced cross section by 44% (34%) at the LHC14 (LHC8), the Tevatron cross section is only 
marginally affected. The relative EW corrections are changed significantly but remain small as 
expected. In contrast, the relative QCD corrections to the total cross sections are reduced from 
~ 40% (30%) at the LHC (Tevatron) to almost zero (10%) after introducing the jet veto proposed 
at the end of Section 14.11 

Predictions for the production cross section with cuts on the transverse momenta or the invari- 
ant mass of the W-pair have been presented already in Section [2j For cross checks and to simplify 
numerical comparisons we show in Tables [3] and H] the cross sections and the corresponding relative 
corrections for different cut values on the W transverse momenta (Table [3]) and on the invariant 
mass of the W-boson pair (Tabled]). While the relative corrections 5 gg due to process (12. 3p are 
negligible at high invariant masses and transverse momenta, we recall that the photon-induced 

3 Although the relative EW corrections hardly depend on the choice of the factorization scale, the QCD cor- 
rections can be significantly stabilized using a scale which is adjusted to the kinematics of the underlying hard 
process. (For a detailed discussion of this important issue, see, e.g., Ref. |40|.) 
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default cuts 


*LO (P b ) 


Sew (%) 


&n(%) 


5 m (%) 


Sqcb (%) 


<5wwv(%) 


LHC8 


23.99 


-0.8 


1.7 


5.8 


1.9 


0.3 


LHC14 


42.39 


-0.9 


1.7 


9.6 


0.5 


0.4 


Tevatron 


7.054 


-0.6 


0.5 


0.5 


11.4 


0.1 



Table 1: Integrated leading-order cross sections and relative corrections for the LHC and the 
Tevatron evaluated with the default setup defined in Section 14.11 



no cuts 


<G (Pb) 


few (%) 


<5 77 (%) 


<W%) 


«S8d (%) 


Swwv(%) 


LHC8 


35.51 


-0.4 


1.5 


5.4 


39.5 


0.3 


LHC14 


75.02 


-0.5 


1.6 


8.1 


44.2 


0.3 


Tevatron 


7.916 


-0.2 


0.5 


0.5 


32.0 


0.1 



Table 2: Total leading-order cross sections and corresponding relative corrections for the LHC 
and the Tevatron without any phase-space cuts. 

channels lead to surprisingly large relative contributions of ~ +20% at Mww ~ 2 TeV, compen- 
sating a significant part of the moderate negative Sudakov-enhanced genuine EW contributions. 
This behaviour may be understood qualitatively recalling that the partonic cross section <r 77 is, 
for sufficiently high energies, given by a constant, <r 77 = 8ira 2 / + O (1/ s) , and strongly peaked 
at small angles, while the corresponding cross section a q q for large s decreases as \n{s/M^)/s. 
Although the NLO QCD corrections turn out to be accidentally small in our default setup, they 
still amount to 20% for large M^\^ and p™*, respectively, exhibiting a similar behaviour at the 
LHC and the Tevatron (see Figs. EJ Hand Tables El SJ) . 



p| ut (GeV) 


<0 (Pb) 


few (%) 


<5 77 (%) 


<W%) 


(%) 


<%wv(%) 


50 


19.80 


-3.0 


2.5 


9.6 


3.5 


0.6 


100 


5.379 


-7.1 


3.8 


9.2 


9.0 


1.0 


250 


35.31 • 10~ 2 


-18.9 


5.8 


7.3 


16.7 


2.3 


500 


23.05 • 10~ 3 


-33.8 


7.0 


3.7 


18.1 


4.0 


750 


33.04 • 10~ 4 


-45.0 


7.9 


2.2 


18.6 


5.3 


1000 


67.04 • 10~ 5 


-53.9 


8.9 


1.4 


18.6 


6.1 


1250 


16.39 • 10~ 5 


-61.4 


10.2 


1.0 


19.5 


6.6 


1500 


44.79 • 10~ 6 


-67.8 


11.8 


0.8 


20.4 


6.9 



Table 3: Integrated leading-order cross sections and relative corrections for different values of 
p™* at the LHC at = 14 TeV, evaluated with the default setup defined in Section |4~T1 
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A^ww (GeV) 


°Z (Pb) 


Sew (%) 


5 77 (%) 


S a (%) 


<^°D (%) 


5wwv(%) 


200 


28.84 


-2.2 


2.2 


8.9 


7.4 


0.5 


300 


9.492 


-4.1 


3.8 


6.4 


14.0 


0.8 


500 


1.841 


-7.5 


7.2 


4.8 


19.2 


1.4 


1000 


12.08 • 10~ 2 


-17.7 


14.4 


2.6 


24.5 


3.2 


1500 


20.37 • 10~ 3 


-25.4 


18.1 


1.4 


23.5 


4.2 


2000 


48.79 • 10" 4 


-31.4 


21.6 


0.9 


23.0 


4.9 


2500 


13.81 • 10" 4 


-36.3 


25.6 


0.6 


22.6 


5.2 


3000 


42.99 • 10~ 5 


-40.6 


30.5 


0.4 


22.4 


5.4 



Table 4: Integrated leading-order cross sections and relative corrections for different values of 
Myfyj at the LHC at yfs = 14TeV, evaluated with the default setup defined in Section [4TT1 

4.3 Transverse-momentum, invariant- mass and rapidity distributions 

In a first step we show the transverse-momentum distributions of W~ and W + (Fig. |5J) and the 
distributions of the invariant mass Mww (Figure [9j) for LHC14, LHC8 and the Tevatron. The 
interpretation of the results is similar to that of the partially integrated cross sections presented in 
Figs. [3] and [4] in Section[2j In particular, we again observe the large Sudakov logarithms resulting 
in corrections of —30% for p^ values around 800 GeV. Again, the correction for the invariant- 
mass distribution is smaller, since large Mww may still involve small momentum transfer t. At 
the LHC, the QCD corrections evaluated with the jet veto proposed in Section 14.11 still reach 
50% (20%) at high p^ (M\vw), while at the Tevatron they do not exceed 20% even at large 
transverse momenta. 

Rapidity distributions of W + and W~ individually are shown in Fig. [TOJ Note that the ra- 
pidity distributions of W + and W - are different at the LHC, as a consequence of the asymmetric 
piece of the differential cross section for uu — > W~W + and dd — > W~W + , and the difference be- 
tween valence- and sea-quark distributions. For the Tevatron one finds d<7w+(y) = dcr w -(— y). 
At the LHC8 (LHC 14), gluon fusion and the 77 process increase the W-pair production rate 
by ~ 5% (10%) and ~ 2% (2%), respectively. The EW corrections are ~ —1%, reflecting the 
small s of typical events and the absence of the large negative Sudakov logarithms. At the Teva- 
tron, only QCD corrections must be considered for this particular observable; EW corrections, 
77 and gg fusion can safely be neglected. 

The rapidity distributions for the LHC 14 and the LHC8 for the subsamples with invariant 
mass Mww > 500 GeV and M\y\y > 1000 GeV are shown in Figs. dU and [I2j respectively. In 
these cases EW corrections and the 77 process are significantly more important. 

Let us finally consider W-pair production at the highest energies accessible at the LHC. 
As already mentioned in Section [21 the distribution in the rapidity difference (for fixed Mww) 
between W + and W~ bosons can be directly related to the angular distribution in the WW rest 
frame and thus can be used to search for anomalous couplings at the TeV scale. In Figs. [T3land[T4l 
we therefore show these distributions for Mww > 1-5 TeV for LHC 14 and Mww > 0.75 TeV 
for LHC8, respectively. At the LHC14, large contributions are visible both from the 77 process 
(+60%) and from the weak corrections (—45%), leading to a strong distortion of this particular 
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Figure 8: Distributions of the W~ transverse momentum at the LHC14 (top), LHC8 (center) 
and the Tevatron (bottom). On the left-hand side, LO contributions due to processes (\2.l$ (qq), 
(|2.2P (ty). and (|2.3j) (gg) are shown. On the right-hand side, corresponding relative corrections 
are presented, normalized to the dominating LO channel (|2.1|) . See text for details. 
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Figure 9: Distributions of the invariant mass of the W-boson pair at the LHC14 (top), LHC8 
(center) and the Tevatron (bottom). On the left-hand side, LO contributions due to pro- 
cesses (|2.ip (gg). (I2.2|) (77). and (I2.3p (gg) are shown. On the right-hand side, corresponding 
relative corrections are presented, normalized to the dominating LO channel (|2.ip . See text for 
details. 
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Figure 10: Distributions of the rapidities of the W-bosons at the LHC14 (top), LHC8 (center) 
and the Tevatron (bottom). On the left-hand side, LO contributions due to processes (j2.1j) (gg). 
(|2.2p (77). and (|2.3|) (gg) are shown. On the right-hand side, corresponding relative corrections 
are presented, normalized to the dominating LO channel (12. ip . See text for details. 
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Figure 11: Distributions of the rapidities of the W-bosons at the LHC 14 (top) and LHC8 (bot- 
tom) for Mww > 500 GeV. On the left-hand side, LO contributions due to processes <\2.1^ (qq), 
(|2.2|) (^). and (12.30 (gg) are shown. On the right-hand side, corresponding relative corrections 
are presented, normalized to the dominating LO channel (|2.ip . See text for details. 
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Figure 12: Distributions of the rapidities of the W-bosons at the LHC14 (top) and LHC8 (bot- 
tom) for M\vw > 1000 GeV. On the left-hand side, LO contributions due to processes (|2,l|) (qq). 
(|2.2p (77). and (|2.3|) (gg) are shown. On the right-hand side, corresponding relative corrections 
are presented, normalized to the dominating LO channel (12, ip . See text for details. 
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Figure 13: Distributions of the W-boson rapidity gap at the LHC14 for Mww > 1500 GeV. 
On the left-hand side, LO contributions due to processes flUD (gg), ([2,2)) (77). and (|2.3|) (gg) are 
shown. On the right-hand side, corresponding relative corrections are presented, normalized to 
the dominating LO channel (I2.ip . See text for details. 
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Figure 14: Distributions of the W-boson rapidity gap at the LHC8 for Mww > 750 GeV. On 
the left-hand side, LO contributions due to processes (|2.ip (qq), (|2.2|) (^r/). and (|2.3|) (gg) are 
shown. On the right-hand side, corresponding relative corrections are presented, normalized to 
the dominating LO channel (I2.ip . See text for details. 
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distribution. As already discussed in Section [2] the additional contributions from gg and 77 
fusion are small in the region accessible to the Tevatron and hence will not be discussed further. 

Finally, we compare our predictions to older results obtained in the high-energy approxima- 
tion. Unfortunately, a tuned comparison of our results with the ones presented in Ref. [56] is 
not possible, since in that paper dedicated event-selection cuts on the leptonic decay products of 
the W bosons are applied. However, comparing Fig. 7 (bottom, left) from Ref. [56] with Fig. [5] 
(bottom, right), we find a reasonable agreement for the relative EW corrections within a few 
percent. Here, we assume that each of the charged leptons on average carries away 50% of the 
momentum of the decaying W, and a strong correlation between Ayww an d Ay^7, which seems 
to be justified in case of strongly-boosted Ws. 

5 Radiation of massive vector bosons 

Finally, we discuss the phenomenological effects of additional massive vector bosons in the final 
state produced in the partonic processes 



Although these channels are parametrically suppressed by one order of the EW coupling, collinear 
or soft massive boson radiation may potentially lead to logarithmically enhanced contributions 
and thus needs further investigation. In the numerical analysis we use a simplified approach to 
conservatively estimate the pollution of W-pair production which can be expected from 3-boson 
final states; we treat the additional boson completely inclusively, whereas our default cuts are 
applied to those opposite-sign W-bosons with highest pr- In Fig. [15] we present the correspond- 
ing numerical results for the WWV (V = W, Z) final state as well as the contribution from 
hard-photon radiation (pt,7 > 15 GeV, |y 7 | < 2.5) relative to quark-induced WW production 
at leading order. (The corresponding numbers can be found in Tables [3] and U] respectively.) 
The relative corrections due to massive boson radiation are below 5% even for large transverse 
momenta and invariant masses and therefore of minor importance. Remarkably enough the con- 
tribution from hard photon radiation with these cuts (which is included in the cross sections 
discussed in the previous chapter) is numerically close to the one from additional massive boson 
radiation. 

6 Conclusions 

We have calculated the full NLO EW corrections to W-boson pair production at hadron colliders 
and present predictions valid in the full energy reach of the LHC. At large parton CM energies, 
the relative corrections are dominated by the well-known universal Sudakov logarithms that lead 
to substantial negative contributions, while the QCD corrections turn out to be moderate after 
application of a dynamic jet veto. 

As a surprising new result we find that photon-induced contributions can be of the same 
size as the genuine EW corrections at high energies and moderate scattering angles and thus 
must not be neglected when predicting W-boson pair production at highest energies. At low 
scattering angles, however, the 77 channel dominates the EW contributions at 0(a 3 ) even at 
moderate energies. In the future, the leptonic decays of the W bosons should be included 



mdj -> W~W + W+ 
Uidj -¥ W~W+W~ 
qq -> W~W+Z. 



(5.1) 
(5.2) 
(5.3) 
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Figure 15: Left: Total LO cross sections for WW, WWV (V = W, Z) and WW7 production; 
Right: Corrections due to WWV and WW7 relative to quark-induced W-pair production. See 
text for details. 

properly, allowing for a realistic event definition to match the increasing accuracy of future LHC 
measurements. 
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